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Abstract 



The behavior of the proton magnetic form factor is studied within the mod- 
ified hard scattering picture, which takes into account gluonic radiative cor- 
rections in terms of transverse separations. We parallel the analysis given 
previously by Li and make apparent a number of serious objections. The 
appropriate cut-off needed to render the form-factor calculation finite is both 
detailed and analyzed by considering different cut-off prescriptions. The use 
of the maximum interquark separation as a common infrared cut-off in the 
Sudakov suppression factor is proposed, since it avoids difficulties with the 
a s -singularities and yields a proton form factor insensitive to the inclusion of 
the soft region which therefore can be confidently attributed to perturbative 
QCD. Results are presented for a variety of proton wave functions including 
also their intrinsic transverse momentum. It turns out that the perturbative 
contribution, although theoretically self-consistent for Q 2 larger than about 

6 GeV 2 to 10 GeV 2 , is too small compared to the data. 
12.38.Bx, 12.38.Cy, 13.40.Gp 
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I. INTRODUCTION 



The proton magnetic form factor at large momentum transfer has been extensively ana- 
lyzed within perturbative Quantum Chromodynamics (pQCD) over the last decade [l -10 



The theoretical basis of these calculations is the hard scattering formula |I;11] in which 
the proton form factor is generically expressed as a convolution of a hard-scattering ampli- 
tude Th and proton distribution amplitudes (DA) $ which represent valence quark Fock 
state wave functions integrated over quark transverse momenta (defined with respect to the 
momentum of their parent proton) : 

G M (Q 2 ) = f\dx\ (\dx'}\f N (ii F )\ 2 $V, Hf)T h (x,x',Q, fi F ), (1.1) 

Jo Jo 

where Q 2 is the invariant momentum transfer squared and [dx] = dx\dx2dx^5{l — 
Xi being the momentum fractions carried by the valence quarks. The renormalization scale 
is denoted by \l and the factorization scale by \xp. The latter scale defines the interface 
between soft physics — absorbed in the wave function — and hard physics, treated explicitly 
within pQCD. The dimensionful constant /jy represents the value of the proton wave function 
at the origin of the configuration space and has to be determined nonperturbatively |§f||§. 
The residual (mainly perturbative) scale dependence of Jn and that of the proton DA is 
controlled by the evolution equation [I]]. 

To lowest order the hard scattering amplitude is calculated as the sum of all Feynman 
diagrams for which the three quark lines are connected pairwise by two gluon propagators. 
This allows the quarks in the initial and final proton to be viewed as moving collinearly up 
to transverse momenta of order ftp. It is then easy to show that Tjj ~ ^ tt '^ , wherein a s 
is the running strong coupling constant in the one-loop approximation. 

The Pauli form factor F2 and hence the electric form factor Ge cannot be calculated 
within the hard scattering picture (HSP), since they require helicity-flip transitions which 
are not possible for (almost) massless quarks. These form factors are dominated by sizeable 
higher twist contributions as we know from experiment [|1^,[K|]. Eq. ( |1 . 1| ) is obtained by 
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taking the + component of the electromagnetic vertex and represents the helicity-conserving 
part of the form factor. 

The choice of the renormalization scale in the calculation of the proton form factor is a 
crucial point. Most authors P,|3|,|7j,P,|T0|1 use a constant a s outside the integrals over fractional 
momenta, with an argument rescaled by the characteristic virtualities for each particular 
model DA. Choosing /i that way and using DAs calculated by means of QCD sum rules — 
distributions whose essential characteristic is a strong asymmetry in phase space — results 
for Gm have been obtained |8| JlO|JT^] that compare fairly well with the data [ 12"lJl5 |. On the 



other hand, the so-called "asymptotic" DA [||] $ as = 120 x\x 2 x^ — into which any DA should 
evolve with Q 2 — > oo — yields a vanishing result for G P M . However, for a renormalization 
scale independent of x, large contributions from higher orders are expected in the endpoint 
region, Xi — > 0. Indeed, for the pion form factor this has been shown explicitly, at least for 



the next-to-leading order [16,17|. Such large higher-order contributions would render the 



leading-order calculation useless. A more appropriate choice of the renormalization scale 



would be, e.g., \Jx2x' 2 Q, since such a scale would eliminate the large logarithms arising 
from the higher-order contributions. Unfortunately, this is achieved at the expense that a s 
becomes singular in the endpoint regions. It has been conjectured |TJ that gluonic radiative 
corrections (Sudakov factors) will suppress that a s -singularity and, therefore, in practical 
applications of the HSP one may handle this difficulty by cutting off a s at a certain value, 
typically chosen in the range 0.5 to 0.7. Another, semi-phenomenological recipe to avoid the 
singularity of a s is to introduce an effective gluon mass |31| which cut-offs the interaction 
at low Q 2 values. 

Besides the extreme sensitivity of the form factors on the utilized DA and besides the 
problem with higher-order contributions and/or the singularity of a s , there is still another — 
perhaps more fundamental — difficulty with such calculations. Indeed, the applicability of 
( |1.1| ) at experimentally accessible momentum transfer, typically a few GeV, is not a priori 
justified. It was argued by Isgur and Llewellyn-Smith [JTS] and also by Radyushkin [2D] that 



the HSP receives its main contributions from the soft endpoint regions, rendering the pertur- 



bative calculation inconsistent. Recently, this criticism has been challenged by Sterman and 



collaborators pl|-|23|j. Based on previous works by Collins, Soper, and Sterman [£4j, they 
have calculated Sudakov corrections to the hard-scattering process taking into account the 
conventionally neglected transverse momentum, k±, of the quarks. The Sudakov corrections 
damp those contributions from the endpoint regions in which transverse momenta of the 
quarks are not large enough to keep the exchanged gluons hard. Moreover, as presumed, the 
Sudakov corrections cancel the a s - singularity without introducing additional ad hoc cut-off 
parameters as for instance a gluon mass. Thus the modified HSP provides a well-defined 
expression for the form factor which takes into account the perturbative contributions in a 
self-consistent way, even for momentum transfers as low as a few GeV. 



However, an important element has not been considered in the analyses of Refs. |f22| , |23 
This concerns the inclusion of the intrinsic transverse momentum of the hadronic wave 
function. As it was recently shown by two of us [25] for the case of the pion form factor, the 



inclusion of the transverse size of the pion extends considerably the self-consistency region 
of the perturbative contribution down to values of momentum transfer unreachable by the 
Sudakov corrections alone. On the other hand, the incorporation of the /^-dependence 
leads to a substantial decrease of the magnitude of the (leading-order) pion form factor. 
Unfortunately, a clear-cut comparison with the available data is not possible because of 
their low quality and the uncertainty in the determination of the pion-nucleon coupling 
constant |26|. Nevertheless, it seems reasonable to expect that the pion form factor receives 



considerable soft contributions in the presently accessible GeV region. 

The aim of the present paper is to perform an analysis for the proton form factor within 
the modified HSP. One of our objectives is to critically examine Li's approach |23j and to 



enlarge the theoretical framework by including the intrinsic fc^-dependence of the proton 
wave function. At the same time we want to clarify several technical points, which are 
absent in the pion case and are first encountered in the more complicated calculation of the 
proton form factor. 

The purpose of our analysis is to investigate how reliably the perturbative contribution 



to the proton form factor can be calculated and to answer the question whether there is 
a proton wave function — modeled on the basis of QCD sum rules j|,|8| — which is capable 
of providing, in a theoretically self-consistent way, a good agreement with the data within 
the modified HSP. It is clear that being able to identify the leading-order perturbative 
contribution reliably allows us to estimate the size of soft contributions to the proton form 
factor, contributions which are not accounted for in the modified HSP formalism. [Note 
that the fcj_-dependent effects taken into account in the modified HSP represent also soft 
contributions of higher-twist type.] 

Sudakov suppression (which can be viewed as the perturbative part of the transverse 
wave function) and intrinsic fc^-dependence of the wave function may also have a lot of 
interesting consequences in other exclusive reactions. Thus, for instance, Sotiropoulos and 
Sterman |27| have applied these elements to near-forward proton-proton elastic scattering 



claiming that their interplay drives the transition of the fixed s differential cross section from 
the t~ s behavior at moderate t to the t~ 10 behavior at larger t, as predicted by dimensional 
counting rules PS . 



The outline of the paper is as follows. In Sec. II we discuss the proton wave function. The 
modified HSP is treated in Sec. III. The discussion of the infrared (IR) cut-off prescription 
in the Sudakov factor and its effect on the a s -singularities is given in Sec. IV. The numerical 
results are presented in Sec. V and our conclusions are contained in Sec. VI. 

II. THE PROTON WAVE FUNCTION 



Similarly to Sotiropoulos and Sterman [p7], we write the valence quark component of the 



proton state with positive helicity in the form 
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\P,+>= f [dx] J [d 2 k ± ] { ^ 123 M a ^" + ^ 213 M a 2. 



0.20.3 



- (*i32 + ^23i)M a + ^_ a:i }e aia2a3 , (2.1) 

where we have assumed the proton to be moving rapidly in the 3-direction. Hence, the 
ratio of transverse to longitudinal momenta of the quarks is small. The measure over the 
transverse momentum integration is defined by 

[^J = 777^1 <^ (2) (£ k ± ) d 2 k ±l d 2 k ±2 d 2 k ±3 . (2.2) 

(167T 3 ) \i=l / 

In the zero binding energy limit, which is characteristic for the parton picture, one has 

xi + x 2 + x 3 = 1 and k ±1 + k± 2 + k± 3 = 0. (2.3) 

The three quark state with helicities Ai, A 2 , A 3 and colors a 1; a 2 , a 3 is given by 

•^aIaoM = . 1 \u ai ;xx, k±i, Ai > \u a2 ; x 2 , k± 2 , A 2 > \d a3 ;x 3 , k ±3 , A 3 > . (2.4) 

3 y/XlX 2 X 3 

Since the orbital angular momentum is assumed to be zero, the proton helicity is the sum 
of the quark helicities. The quark states are normalized as follows: 

< ?«;; 4 k'^Klqa^x^ k ±i X t > = 2x l (2^f5 <a 5 Kx 5(x' i - x { )5(k' ±i - k ±i ). (2.5) 

From the permutation symmetry between the two u quarks and from the requirement that 
the three quarks have to be coupled to give an isospin 1/2 state it follows that Eq. fl2.1|) can 
be expressed in terms of only one independent scalar function 0. In the sequel, \1/ denotes 
the momentum space wave function. 

The subscripts on ^ refer to the order of momentum arguments, for example 
^123(2, k±) = ^(xi, k±i, x 2 , k± 2 ; x 3 , k± 3 ). Note that, in general, the wave function depends 
on the factorization scale /xp. We make the following convenient ansatz for the wave function: 

*i23 0e,£l) = —== f N {^ F )^{x,ii F )n{x,k ± ). (2.6) 

The distribution amplitude $ = V — A (in the notation of H) is defined in such a way that 



1 

[dx]$x2a(x,fi F ) = l, (2.7) 

o 

where an obvious abbreviated notation has been introduced. The DA can be expressed in 
terms of the eigenfunctions of the evolution equation [|IJ, $ ra (xj), which are linear combina- 
tions of Appell polynomials. Then the proton DA can be cast into the form 

fciasOr, a*) = *5s(s) £#n [^^Y^ ( 2 - 8 ) 

where the notations of |7j are adopted. $"23 ( x ) * s the asymptotic DA mentioned in the 
Introduction. The exponents j n , driving the evolution behavior of the DA, are related to the 
anomalous dimensions of trilinear quark operators with isospin 1/2 (see |[29|| ) and resemble 
the b n in the Brodsky-Lepage notation Because they are positive fractional numbers 
increasing with n, higher-order terms in ( |2.8|) are gradually suppressed. The constants j n 
are given in Table 1; /3o = 11 — 2rif/3 = 9 for three flavors. 

Constraints on the DA are obtained implicitly by restricting their few first moments 
within intervals determined from QCD sum rules [@,|]|§, which are evaluated at some self- 
consistently determined normalization point fip of order 1 GeV (see, e.g., 0): 

K ( " 1W) W= [\dx]x^x^x^^ 123 (x,^) (2.9) 







In most model calculations, mentioned above, the moment constraints provided by QCD 
sum rules are used to determine the first five expansion coefficients B n , where B = 1 due to 
normalization (2.7). However, since the moments are burdened by errors, these expansion 
coefficients — although mathematically uniquely determined by the moments of correspond- 



ing order [T^ — in practice their numerical values cannot be fixed precisely giving rise to 
different options for the proton DA. In our calculation of form factors we employ amplitudes 
complying with the Chernyak-Ogloblin-Zhitnitsky (COZ) sum-rule moment constraints. It 



was shown in [ |30| , |3T| that such amplitudes constitute a finite orbit in the (-B4, R = {G^l/Gfj^) 



plane ranging from COZ-like amplitudes M with R < 0.5 to the recently proposed [TIJ het- 
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erotic one with R m 0.1. For the convenience of the reader, the QCD sum-rules constraints 
and the expansion coefficients B n of selected model amplitudes are compiled in Table 1. 

The fcj_-dependence of the wave function is contained in the function Q which is normal- 
ized according to 

"[d 2 k ± }n 123 (x,k ± ) = l. (2.10) 



Due to Q2.7D and fl2.10|) , /n is the value of the DA at the origin of the configuration space. 



Its evolution behavior is given by 

and its value has been determined to be /at(^o) = (5-0 ± 0.3) x 10 _3 GeV 2 @||. 

In Eq. (|2.6|) \I/ represents the soft part of the proton wave function, which results by 
removing the perturbative part and absorbing it into the hard-scattering amplitude T#. 
The perturbative tail of the full wave function behaves as l/k\ for large k±, whereas the 
soft part vanishes as l/k\ or faster. The nonperturbative or intrinsic fc^-dependence of 
the soft wave function, being related to confinement, is parametrized as a simple Gaussian 
according to 

3 



^123 fa, k±) = (167r ) exp 

XiX 2 X 3 



-a 2 ^2 k±i/ x i 

i=l 



(2.12) 



This parametrization of the intrinsic /^-dependence of the wave function, which is due to 
Brodsky, Huang, Lepage, and Mackenzie p2| , seems to be more favorable than the standard 



form of factorizing x- and /c^-dependencies. At least for the case of the pion wave function, 
this has recently been effected by Zhitnitsky |33| on the basis of QCD sum rules. He finds 



that a factorizing wave function is in conflict with some general theoretical constraints which 
any reasonable wave function should comply. Zhitnitsky's QCD sum-rule analysis of the 
pion wave function seems to indicate that the fc^-distribution may also show a double-hump 
structure, which means that small and large values of k± are favored relative to intermediate 
values. It is likely that the proton wave function may exhibit a similar behavior, though 
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this kind of analysis has yet to be done. For the purposes of the present work we ignore this 
possibility. 

In ( p.!2|) the parameter a controls the root mean square transverse momentum (r.m.s.), 
(k^) 1 / 2 , and the r.m.s. transverse radius of the proton valence Fock state. From the known 
charge radius of the proton, we expect the r.m.s. transverse momentum to be larger than 
about 250 MeV. The actual value of (Ar[) 1//2 may be much larger than 250 MeV, e.g., 600 MeV 
or so. Indeed, Sotiropoulos and Sterman p7| show that application of the modified HSP 



to proton-proton elastic scattering leads to an approximate t _8 -behavior of the differential 
cross section at moderate \t\. The behavior da/dt ~ t~ 10 , predicted by dimensional counting, 
appears only at very large \t\. At precisely which value of \t\ the transition from the t~ 8 
to the t~ 10 behavior occurs, depends on the transverse size of the valence Fock state of the 
proton. Since the ISR [34]] and the FNAL |35j data are rather compatible with a t _8 -behavior 



of the differential cross section, Sotiropoulos and Sterman conclude that the transverse size 
of the proton is small, perhaps < 0.3 fm. Correspondingly, the r.m.s. transverse momentum 
is larger than 600 MeV. It is worth noting that such a value is supported by the findings of 



the EMC group [3~6| in a study of the transverse momentum distribution in semi-inclusive 
deep inelastic fip scattering. A phenomenologically successful approach to the HSP, in which 
baryons are viewed as bound states of a quark and an effective diquark, also uses a value 
of this size for (k 2 ^ 1 ^ 2 |37|-p9|. There is a second constraint on the wave function, viz. the 



probability for finding three valence quarks in the proton: 



vra) 4 / [dx] K 1MK " *^">*^"> < L ( 2 .13) 



3 Jo X1X2X3 

In our numerical analysis to be presented in Sec. 5, we make use of two different values of 
the r.m.s. transverse momentum, namely, one which is obtained by the requirement P 3q = 1 
for a given wave function. [This corresponds to the minimum value of the r.m.s. transverse 
momentum.] The other option for the r.m.s. transverse momentum we consider is the rather 
large value of 600 MeV. In the latter case, the probability for the valence quark Fock state 
depends on the wave function. 
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III. THE MODIFIED HARD SCATTERING PICTURE 



Following Li [2^, we write the proton form factor in the form 

G M (Q 2 ) = — J o [dx}[dx'} J [GpfcJ[d 2 /4] ^T^x,^^,^ (3.1) 

Note, however, that our notation is slightly different compared to that of Li. Making use of 
the symmetry properties of the proton wave function under permutation, the contributions 
from the 42 diagrams involved in the calculation of the proton form factor in lowest order 
can be arranged into two reduced hard scattering amplitudes of the form 



■Hi 



■a 



F T 



(47ra s (/j)) : 



;i - xi)(i - x'oq 2 + (k ±1 - k' ±1 y x 2 x' 2 q* + (k ±2 - k' ±2 y 



(3.2) 



Th 2 -;C f 



(47ra,(/i)) s 



X!x[Q 2 + (k±i - k' ±1 ) 2 ] [x 2 x' 2 Q 2 + (k ±2 - k' ±2 ) 



(3.3) 



where Cp = 4/3 is the Casimir operator of the fundamental representation of SU(3) C . In 
the hard scattering amplitudes only the fc^-dependence of the gluon propagators is included, 
whereas that of the quark propagators has been neglected. It is expected that this technical 
simplification introduces only a minor error of about 10% in the final result. For the case of 
the pion form factor this has been explicitly demonstrated by Li [23 . 

The functions Yj in ( [3.1|) are short-hand notations for linear combinations of products of 
the initial and final state wave functions ^ijk^*!ji k >, weighted by x^-dependent factors 
from the fermion propagators, namely: 



; arising 



Y 1 = 



4^23*123 + 4^ 32 * 13 2 + 



231^231 + ^321^321 



+ 2^31^132 + 2^32^231 + 2tf£' 21 tf i 23 + 2#f 23 # 321 } 



(3.4) 



Yo 



2(l-x 2 )(l-x' 1 
1 



- {3^32^132 - ^231*231 ~ *231*132 ~ ^32*231 } 



(1 - x 3 )(l - x[) 



j- {4^ 21 ^ 321 + tff 23 * 123 + 2** , 21 * 123 + 2^' 2 3^3 21 } 



(3.5) 
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Ignoring the transverse momenta in the hard scattering amplitudes ( |3.2|) and ( |3.3|) , and 
inserting ( |2.6| ) and ( |2.10| ), one arrives at the standard HSP result for the magnetic form 



factor. Although this expression is correct in the asymptotic momentum domain, the trans- 
verse degrees of freedom are an essential ingredient of the formalism and neglecting them 
leads to inconsistencies in the endpoint regions, where one of the fractional momenta Xj 
or x\ tends to zero. After all, it is precisely this approximation that is responsible for the 
inconsistencies mentioned in the Introduction. The power of combining the transverse mo- 
mentum dependence of the hard scattering amplitude and radiative corrections in the form 



of Sudakov form factors was realized by Sterman and collaborators |2l|-p3|. Ultimately, 
it leads to a suppression of contributions from the dangerous soft regions, where both the 
longitudinal and transverse momenta of the quarks are small. 

In order to include the Sudakov corrections, it is advantageous to reexpress Eq. ( |3.1|) in 
terms of the variables bi, which are canonically conjugate to k±i and span the transverse 
configuration space. Then 

Gm(Q 2 ) = yJ l dx W\ J (4^)2 (4^)2 E f i(x, b, Q, n)Yj{x, x', b, » F ) e^, (3.6) 

where the Fourier transform of a function f(k±) = f(k±i, k± 2 ) is defined by 

f(b) = j^y J d 2 k ± id 2 k ±2 exp{-ib v k ±1 - ib 2 - k± 2 }f(k). (3.7) 

Since the hard scattering amplitudes depend only on the differences of initial and final 
state transverse momenta, there are only two independent Fourier-conjugate vectors b\ (= 
b[) and b 2 (= b' 2 ). They are, respectively, the transverse separation vectors between quarks 
1 and 3 and between quarks 2 and 3. Accordingly, the transverse separation of quark 1 and 
quark 2 is given by 

h = b 2 -h. (3.8) 

[Note that Sotiropoulos and Sterman J^j] define the transverse separations in a cyclic way 
which results in the interchange b\ < — > —b 2 , as compared to our definition.] 
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The fact that there are only two independent transverse separation vectors is a conse- 
quence of the approximation made in the treatment of the hard scattering amplitudes ( |3.2| ) 
and ( |3.3| ) which disregards the fcj_-dependence of the quark propagators. This approxima- 
tion is justified by the enormous technical simplification it entails, given that the thereby 
introduced errors are very small. Then by virtue of rotational invariance of the system with 
respect to the longitudinal axis, the form factor (|3.6| ) can be expressed in terms of a seven- 
dimensional integral instead of an eleven- dimensional one. Physically, the relations b\ = b[, 
b 2 = b 2 mean that the physical probe (i.e., the photon) mediates only such transitions from 
the initial to the final proton state, which have the same transverse configurations of the 
quarks. 

The Fourier-transformed hard scattering amplitudes appearing in Eq. ( |3.6|) read 

Ti = I C F a s {t n )a s {t 12 )K - rr a )(l - x' x )Qb^j K Ux~^ 2 Qb 2 \ , (3.9) 



T 2 = ~C F a s (t 21 )a s (t 22 )K [yxx^QbA K \Jx 2 x 2 Qb 2 j , (3.10) 

where Kq is the modified Bessel function of order and bi denotes the length of the corre- 
sponding vector. We have now chosen the renormalization scale in such a way that each hard 
gluon carries its own individual momentum scale tji as the argument of the corresponding 
a s . The tji are defined as the maximum scale of either the longitudinal momentum or the 
inverse transverse separation, associated with each of the gluons: 



yJ(l-X 1 )(l-X , 1 )Q,l/b 1 



tn= max 

t 2 i= max \Jxxx' x Q, l/b x 
tu= t 2 2 = max 



x 2 x' 2 Q, l/b 2 



(3.11) 



One may think of other choices. However, they are not expected to lead to very different 



predictions for the form factor [23 



The quantities Yj contain the same combinations of initial and final state wave func- 
tions as those in ( |3.4|) and ( p.5|) , the only difference being that now the products 
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ty*!j, k ,^ijk are replaced by corresponding products of Fourier-transformed wave functions: 
^f*!j, k ,(x' ,b, fj,F)^fijk(x,b, Hf)- Using (|2.6|) and ( |2.12| ), the Fourier transform of the wave 
function reads 



^I2s(x, b, /i F ) 



/n (Vf) $123 0, /U F )^i23 {X, b) , 



(3.12) 



where the Fourier-transform of the k±- dependent part is given by 

1 



^i2s(x,b) = (Air) exp 



4a 2 



xix 3 b 1 + x 2 x 3 b 2 + xix 2 b 3 



(3.13) 



The exponentials e •» in ( |3.6| ) are the Sudakov factors, which incorporate the effects of 
gluonic radiative corrections. Because of this, ( |3.6| ) is not simply the Fourier transform of 
( 3.1|) but an expression comprising an additional physical input. Thus (|3.6|) may be termed 
the "modified hard-scattering formula". On the ground of previous works by Collins and 
Soper |24]| , Botts and Sterman |21] have calculated a Sudakov factor using resummation 
techniques and having recourse to the renormalization group. They find Sudakov exponents 
of the form 



1=1 L 

3 

+ E 

i=i 



dfx 



s{x h b h Q)+ I _ —i q (g(v 2 )) 

l/bi /i 



s(x' l ,b l ,Q)+ f 3 !%M?)) 

Jl/bi /! 



(3.14) 



wherein the Sudakov functions s(£j, bi, Q) are given by 

ste , 5l , Q)= ^ 41n (lU^f*-i 



2A 



\bi 



AM fa 
16/3? 



ln(26,) + 1 ln(2ft) + 1 



ft 



4(2) 4(1) 

^_ _ ^_i n fe 2 'r- 1 /2 



In 



32/3? 



In 2 (2g0-ln 2 (2^ 



(3.15) 



Here ^ = X\ or a;J (/ = 1, 2, 3) and the variables q and 6 are defined as follows: 



q t = ln[&Q/(V2A QCD )] 



(3.16) 
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k = ln[l/6,A QCD ]. (3.17) 



The coefficients A® and /3, are 



Ad) 4 a(2) 67 1 2 10 8 _ . /I _ 

= -, = 7T 2 n,/ + -ftln -e 7 

3 9 3 27 ' 3 V2 

A = i*^, A = Hz^5£, ,3,8) 

where to/ is the number of quark flavors and 7 = 0.5772 ... is the Euler-Mascheroni constant. 
In the sequel n/ = 3 is used. 7 9 = — — + O(a^) is the anomalous quark dimension in the 
axial gauge @ . 



The Sudakov function, s(&, bi,Q), in ( p.!5|) takes into account leading and next-to- leading 



gluonic radiative corrections of the form shown in Fig. |l[ The quantities bi (I = 1,2,3) 
are infrared cut-off parameters, naturally related to, but not uniquely determined by the 
mutual separations of the three quarks f24|| . A physical perspective on the choice of the 



IR cut-off is provided by the following analogy to ordinary QED. One expects that because 
of the color neutrality of a hadron, its quark distribution cannot be resolved by gluons 
with a wave length much larger than a characteristic quark separation scale; meaning that 
long wave length gluons probe the color singlet proton and hence radiation is damped. 
Radiative corrections with wave lengths between the IR cut-off and an upper limit (related 
to the physical momentum Q) yield to suppression; it is understood that still softer gluonic 
corrections are already taken care of in the hadron wave function, whereas harder gluons 
are considered as part of T# . 

Different choices of the IR cut-off have been used in the literature: Thus, Li |23| chooses 



bi = b\ (this choice hereafter is termed the "L" prescription), whereas Hyer j|T| in his 



analysis of the proton-antiproton annihilation into two photons and of the time-like proton 
form factor as well as Sotiropoulos and Sterman |?7| take b\ = b 2 , 62 — b\, b 3 = 6 3 (this 
choice is denoted the "H-SS" prescription). Still another possibility, and the one proposed 
in the present work for reasons that will be explained below is to use as IR cut-off the 
maximum of the three interquark separations, i.e., to set 
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b = max{6i, b 2 , b 3 } = b x = b 2 = b 3 . (3.19) 

This choice, designated by "MAX" , is analogous to that in the meson case, wherein the 
quark- ant iquark distance naturally provides a secure IR cut-off. The specific features of each 



particular cut-off choice will be discussed in detail in Sec. IV. 



The integrals in ( |3.14| ) arise from the application of the renormalization group equation 



(RGE). The evolution from one scale value to another is governed by the anomalous dimen- 
sions of the involved operators. The integrals combine the effects of the application of the 
RGE on the wave functions and on the hard scattering amplitude. The range of validity 



of (|3.15|) for the Sudakov functions is limited to not too small bi values. Whenever 1/bi 
is large relative to the hard (gluon) scale £iQ, the gluonic corrections are to be considered 
as higher-order corrections to Tjj and hence are not contained in the Sudakov factor but 



are absorbed in T#. For that reason, Li [23[ sets any Sudakov function s(£i,bi,Q) equal to 
zero whenever ^ < y/2 /(Qbi). Moreover, Li holds the Sudakov factor e~ Sj equal to unity 
whenever it exceeds this value, which is the case in the small ^-region. Actually, the full 



expression ( |3.14j ) shows in this region a small enhancement resulting from the interplay of 



the next-to-leading logarithmic contributions to the Sudakov exponents and the integrals 
over the anomalous dimensions. We follow the same lines of argument in our analysis. 

The IR cut-offs 1/6; in the Sudakov exponents mark the interface between the nonpertur- 
batively soft momenta, which are implicitly accounted for in the proton wave function, and 
the contributions from soft gluons, incorporated in a perturbative way in the Sudakov fac- 
tors. Obviously, the IR cut-off serves at the same time as the gliding factorization scale /xp to 
be used in the evolution of the wave function. For that reason, Li |23| as well as Sotiropoulos 
and Sterman [^7]] take fip = min{l/fy}. The "MAX" prescription (|3.19|) , adopted in the 
present work, naturally complies with the choice of the evolution scale proposed in |23"1 , |27 . 
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IV. DISCUSSION OF THE ^-SINGULARITIES 



It is well known that the inclusion of an x-dependent renormalization scale in the argu- 
ment of a s within the standard HSP of Brodsky-Lepage JI] presents the difficulty that the 
value of a s becomes singular in the endpoint regions. To render the form factors (Eq. (|1 . 1| ) ) 



finite, additional external parameters, like an effective gluon mass [I8| or a cut-off prescrip- 
tion have to be introduced. Technically, such parameters play the role of IR regulators 
serving to regularize one of the gluon propagators, which may become soft along the bound- 
aries of phase space (see, e.g., [[7]). One of the crucial advantages of the modified HSP, 
proposed by Sterman and collaborators [pl~H 23H, is that there is no need for external regula- 
tors because the Sudakov factor may suppress the singularities of the "bare" (one-loop) a s 



inherently. Indeed in the pion case, it was shown |22j that the transverse quark- ant iquark 
separation is tantamount to an IR regulator which suffices to cancel all singularities from 
the soft region. 

Concerning the proton form factor, the situation is much more complicated because more 
scales are involved and hence the choice of the appropriate IR cut-off parameters b\ is not 
obvious, as discussed in Sec. III. As we shall effect in the following, the cancellation of the 
a s -singularities by the Sudakov factor depends sensitively on that particular choice. 

In Fig. we display the exponential of the Sudakov function exp[— s(£i, bi, Q)] for 
Q = 30 Aqcd by imposing Li's requirement |23|| : s(£i,bi,Q) = whenever £j < y/2/Qbi. 
Ultimately, the cancellation of the (^-singularities relies on the fact that whenever one of 
the a s tends to infinity (owing to the limit tji — > Aqcd), the Sudakov factor e - ^ rapidly 
decreases to zero. As it can be observed from Fig. [| this is not the case in the region deter- 
mined by & < v^Aqcd/Q and simultaneously 6;Aqcd — ► 1, where exp[— s(£i, bi, Q)) is fixed 
to unity. In the pion case this does not matter, since the other exp[— s(l — £,b,Q)] — > 
faster than any power of 1u[1/(6Aqcd)] and, consequently, the Sudakov factor drops to zero. 
In contrast, the treatment of the proton form factor is more subtle. In that case, e _Sj does 
not necessarily vanish fast enough to guarantee the cancellation of the a s -singularities. This 
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can be illustrated by the following configuration: if, say, x\ < v^Aqcd/Q and &iAqcd —* 1 
then X2 + X3 ~ 1 and X2 can have any value between and 1 — x\. Since &2Aqcd is unre- 
stricted within the limits and 1, the corresponding exponentials of the Sudakov functions 
exp[— s(x2, 62, Q)} and exp[— s(x3, 63, Q)} do not automatically fall off to zero in order to yield 
sufficient suppression of the a s -singularities, unless all three hi are coerced to be equal. If 
the three hi are allowed to be different, then the Sudakov factor provides suppression only 
through the contributions of the anomalous dimensions. According to the "L" and "H-SS" 
prescriptions, which, in general, allow for different 6/ in the Sudakov functions, the integrand 
in (|3.6|) has singularities behaving as 

~ln(^ — \ (4.1) 

V^Aqcd/ 

for 6jAqcd — 1 and xi hold fixed. The maximum degree of divergence is given by 

K = J it + 2 ^ max " 2 ) + (4 ' 2) 
where the first term 4/3 comes from the evolution of f^, fl2.11| ) and the constant 7 max is 



related to the anomalous dimension driving the evolution behavior of the proton DA, see 
( p.8| ) and Table 1: 7 max is the maximum value of the {7 n } within a given polynomial order of 
the expansion of the DA. We reiterate that the 7 n are positive fractional numbers increasing 
with n. Thus the singular behavior of the integrand becomes worse as the expansion in 
terms of Appell polynomials extends to higher and higher orders. The term —2 in ( |1.2|) 
stems from the integrations over the anomalous dimensions in the Sudakov factor e _5j (see 
( |3.14j )). Finally, the term 1 originates from that a s (tjk) which becomes singular in ( |3.6|) , c.f., 
( [3. 11| ). Which one of the a s couplings becomes actually singular, depends on the prescription 
imposed on the IR cut-off parameters 6;. The integral ( |3.6| ) does not exist if 7 max > |. As 
Table 1 reveals, this happens already for proton DAs which include Appell polynomials of 
order 1, i.e., for all DAs except for the asymptotic one: $ as = 120x1X2X3. Thus application 
of the "L" and "H-SS" prescriptions on the choice of the IR cut-off parameters hi to the 
proton form factor entails the modified HSP to be invalid. In view of these results, Li's 



analysis of the proton form factor [^| seems to be seriously flawed 



A simple recipe to bypass the singular behavior of the integrand is to ignore completely 



suggested to take for the factorization scale fip — max (1/6;) . In this case, the 7 max appears in 
( |4.2| ) only if all three bi tend to 1/Aqcd at once. But then at least one of the exp[— s(£j, bi, Q)] 



ization scale avoids singularities of the form ([O]), but seems to us physically implausible. 
Since he only presents numerical results for the proton form factor in the time-like region, 
we cannot compare with his results directly. 

Another option, and actually the one proposed in this work, is to use a common IR 
cut-off not only for the evolution of the wave function but also in the Sudakov exponent. 
For a common cut-off b, the Sudakov factors always cancel the a s -singularities; if, for a given 
I, we are in the dangerous region, £j < v^Aqcd/Q, 6Aqcd — * 1> at least one of the other 
two Sudakov functions lies in the region & > a/2Aqcd/Q, 6Aqcd — ► 1 (I' 7^ I) and therefore 
provides sufficient suppression, as outlined above. In particular, we favor b = max{6/} as 
the optimum choice ("MAX" prescription), since it does not only lead to a regular integral 
but also to a non-singular integrand. The Sudakov factor e _Sl subject to the "L" and 
"MAX" prescriptions is plotted for a specific quark configuration in Fig. [| This figure 
makes it apparent that the Sudakov factor in connection with the "MAX" prescription 
is unencumbered by singularities in the dangerous soft regions. As a consequence of the 
regularizing power of the "MAX" prescription, the perturbative contribution to the proton 
form factor ( |3.6| ) saturates in the sense that the results become insensitive to the inclusion 
of the soft regions. A saturation as strong as possible is a prerequisite for the self-consistency 
of the modified HSP, as will be discussed in Sec. V. 

To demonstrate the amount of saturation, we calculate the proton form factor through 
( P-6| ), employing a cut-off procedure to the bi- integrations at a maximum value b c . In Fig. |] 
the dependence of Gm on b c for the three choices, labeled: "L", "H-SS", and "MAX" 
is shown using, for reasons of comparison with previous works, the COZ DA and ignoring 
evolution. [Evolution has been dispensed with to avoid the concomitant singularity in Q 4 Gm 



the evolution of the DA or to "freeze" it at any (arbitrary) value larger than Aqcd- Hyer f41| 




Apparently, Hyer's choice of the factor- 
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as 6 c Aqcd — » 1 when imposing the "L" and "H-SS" prescriptions.] As one sees from the 
figure, the "MAX" prescription leads indeed to saturation; the soft region 6 c Aq CD > 0.7 
does not contribute to the form factor substantially. In fact, already 50 % of the result are 
obtained from the regions with 6 c Aqcd < 0.48. Note that a s increases to a value of 0.95 
at 6 c Aqcd ~ 0.48. This indicates that a sizeable fraction of the contributions to the form 
factor is accumulated in the perturbative region. 

Unfortunately, this saturation is achieved at the expense of a rather strong damping of 
the perturbative contribution to the proton form factor. Using the two other prescriptions 
("L" and "H-SS") and ignoring evolution, we have found larger results for Gm, but no 
indication for saturation: The additional contributions to the form factor gained this way 
are accumulated exclusively in the soft regions, i.e., for values of 6 c Aq C d near 1. These 
findings are in evident contradiction to Li's results (figure 5 in |23|]) for which an acceptable 
saturation has been claimed. On the other hand, we can qualitatively confirm the saturation 
behavior of the proton form factor calculated by Hyer |4l| in the time-like region. Since we 
regard a saturation behavior as a stringent test for the self-consistent applicability of pQCD, 
calculations which accumulate large contributions from soft regions (large b c ) cannot be 
considered as theoretically legitimate. 

The role of the evolution effect subject to the "MAX" prescription is also exhibited in 
Fig. |J. It shows that the effect of evolution is large, although finite, owing to the strong 
suppression provided by the Sudakov factor. Note that according to our discussion in Sec. Ill, 
the factorization scale is /ip = 1/b. The significant feature of the evolution effect is that 
it tends to neutralize the influence of the IR cut-off. Thus one obtains larger values of the 
proton form factor at the expense of a slightly worse saturation. 

V. NUMERICAL ANALYSIS 

In this section we give numerical results for the proton form factor. In these cal- 
culations we throughout employ the "MAX" prescription with evolution included, using 
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Aqcd =180 MeV and /io = 1 GeV. Before proceeding with the presentation of our final 
results, let us investigate the effect of including the intrinsic transverse momentum in our 
calculations. The fc^-dependence of the proton wave function effectively introduces a con- 
finement scale in the formalism, the importance of which may be appreciated by looking 
at Fig. []. This figure shows results, obtained for the COZ DA without fc^-dependence and 
for two different values of {k\) 1 / 2 . To describe the intrinsic A;j_-dependence, one can use 
( 2.12Q or, after transforming to the transverse configuration space, (|3l3|) . Notice that in 



Li's approach the Gaussian in (|3.13|) has been replaced by unity. The oscillator parameter a 
is determined in such a way that either the normalization of the wave function P^ q is unity 
(resulting into (k\) 1 / 2 = 271 MeV for the COZ DA), or by inputing the value of the r.m.s. 
transverse momentum. In the second case, we use a value of 600 MeV (see the discussion 
in Sec. II), which implies P% q = 0.042. As can be seen from this figure, the predictions for 
the form factor are quite different for the three cases. The intrinsic /cj_-dependence of the 
wave function leads to further suppression of the perturbative contribution, which becomes 
substantial if the r.m.s. transverse momentum is large. On the other hand, this suppres- 
sion is accompanied by an increasing amount of saturation, since also the Gaussian (|3.13p 



suppresses predominantly contributions from the soft regions, viz., the large 6-regions. In 
contrast to the Sudakov factor, however, this suppression is Q-independent. The interplay of 
the two effects, Sudakov suppression and intrinsic transverse momentum, leads to a different 
Q-behavior of the form factor depending on the value of the r.m.s transverse momentum, as 
can be seen from Fig. |5|. The Q-dependence beyond 10 GeV 2 is rather weak, being approx- 
imately compatible with dimensional counting (modulo logarithmic corrections). For very 
large values of Q beyond 1000 GeV 2 the three curves have approached each other within 
10% accuracy. This happens when the Sudakov factor dominates the Gaussian ( |3.13[ ) and 



selects those configurations with small interquark separations. In this region, which one may 
consider as the pure perturbative region, the results for the form factor are independent of 
the confinement scale introduced by the r.m.s. transverse momentum. 

The penalty of the additional suppression of the perturbative contribution caused by the 
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Gaussian ( |3.13|) is mitigated by the advantage that the perturbative contribution becomes 
more self-consistent than by the Sudakov factor alone. This is indicated in the enhanced 
amount of saturation with increasing r.m.s. transverse momentum. Adapting the criterion 



of self-consistency, originally suggested by Li and Sterman |E| for the pion case, namely 
that 50% of the results are accumulated at moderate values of the coupling constant, say, 
a 2 < 0.5, we find self-consistency for Q 2 « 7 GeV 2 (for the COZ DA). 

Finally, in Fig. |6|, we demonstrate the effect of different proton DAs on the form factor. 



To this end, we investigate a set of 45 DAs [pOipil , which all respect the QCD sum-rule 



constraints |§. The results for the various DAs — or more precisely wave functions, since 
we include their intrinsic transverse momentum dependence — obtained under the "MAX" 
prescription with evolution included, form the shaded area shown in the figure. All wave 
functions are normalized to unity and the corresponding r.m.s. transverse momenta vary 
between 267 MeV and 317 MeV (see Table 1). The theoretical form-factor predictions span 



a "band" congruent to the "orbit" of solutions found in ||3Q|J3T|] . The upper bound of the 
"band" corresponds to the DA COZ up , which yields the maximum value of the form-factor 
ratio {G^l/G^j = 0.4881 in the standard HSP. The lower limit of the "band" is obtained 
using the DA "low" (sample 8 in 0) with \G n M \/G p M = 0.175. Explicitly shown are the 
results for the COZ DA ||, its optimized version (with respect to the sum-rule constraints) 



and the "heterotic" DA, recently proposed by two of us in [10]. We note that the differences 
among these curves practically disappear already at about Q 2 = 80 GeV 2 , despite the fact 



that these amplitudes have distinct geometrical characteristics [|14 . 

Since the true valence Fock state probability is likely much smaller, or invariably the 
r.m.s. transverse momentum larger than of order of 300 MeV, the "band" describes rather 
maximal expectations for the (leading-order) perturbative contributions to the form factor; 
at least for proton wave functions of the type we utilize. Comparison with the experimental 
data reveals that the theoretical predictions amount, at best, to approximately 50% of the 
measured values. This is the benchmark against which we have to discern novelties and 
aberrations. Closing this discussion we note that, since we are calculating only the helicity- 
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conserving part of the current matrix element it is not obvious whether we should compare 
the theoretical predictions with the data for the Sachs form factor Gm or the Dirac form 
factor F\. Therefore we have exhibited in Fig. [] both sets of data [[HJ for comparison. Since 



the two sets of data differ by only 10%, our conclusions concerning the smallness of the 
theoretical results remain unaffected. 

The various model wave functions led to self-consistency of the perturbative contribution, 
i.e., 50% of the results are accumulated in regions where a 2 s < 0.5, in the range of Q 2 between 
6 and 10 GeV 2 . 

VI. SUMMARY AND CONCLUSIONS 

The objective of the present work has been to derive the proton magnetic form factor 
within the modified version (Sec. Ill) of the standard Brodsky-Lepage HSP 0, a scheme 
which takes into account gluonic radiative corrections |3| in terms of transverse separa- 
tions. This is done by incorporating in the formalism the Sudakov factor, calculated by 
Botts and Sterman There are already some interesting applications of the modified 



HSP [^, P^ , Pq , P7| , ^l| -^3| . The significant element of this type of analyses is that the a s - 
singularities, arising from hard-gluon exchange and evolution, can be cancelled without 
introducing free external parameters. We emphasize that in contrast to pure phenomeno- 
logical recipies (e.g., the introduction of a gluon mass), the modified HSP provides an explicit 
scheme how the IR protection of the "bare" a s proceeds through gluonic radiation accumu- 
lated in the Sudakov factor. Thus, in the modified HSP, one may conceive of the (finite) 
IR-protected a s as being the effective coupling. By this procedure the potentially dangerous 
soft regions of momenta are suppressed entailing also a reduction of the perturbative contri- 
bution to the form factor. While in the pion case p2[ , it is fortunate that the cancellation of 



the a s -singularities comes out naturally, Li's approach to the proton form factor leads to 
a lack of complete cancellation of the a s -singularities (see Sec. IV). Without evolution of the 
proton wave function the emerging singularities in (|3.(j|) are still integrable, but logarithmic 
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corrections due to evolution yield ultimately to uncompensated singularities. 

On the grounds of our discussion, we are reasonably confident that Li's treatment can 
be cured within the modified HSP. We suggest to use a a common IR cut-off in the Sudakov 
exponents ( |3.f4| ) and Sudakov functions ( |3.f 5| ): viz., the maximum transverse separation. 



This "MAX" prescription provides sufficient IR protection, since even with evolution, the 
integrand in (|3.6| ) remains finite. A significant feature of this treatment is that the proton 
form factor saturates, i.e., it becomes insensitive to the contributions from large transverse 
separations. The other choices of the IR cut-off ("L", "H-SS"), we have discussed, do not 
lead to saturation. 

However, this reliable saturation and IR protection of the form factor is achieved at 
the expense of a strong reduction of the perturbative contribution to the form factor. The 
damping of the proton form factor becomes even stronger if one takes into account the 
intrinsic transverse momentum dependence of the proton wave function (see Fig. f| and 
Fig. |5|). This has been done by assuming a non-factorizing x and fc^-dependence of the 
wave function of the Brodsky-Lepage-Huang-Mackenzie E| type and fixing the value of 



(fc^.) 1 / 2 either via the valence quark probability P^ q or by inputing the value 600 MeV by 
hand ||27|| . A remarkable finding is that the form factor calculated within the modified HSP, 
appropriately extended to include the intrinsic transverse momentum of the proton wave 
function, shows only a mild dependence on the particular model DA. 

The perturbative contribution to the form factor becomes self-consistent in all cases 
for momentum transfers larger than 6 to 10 GeV 2 . The actual value of the onset of self- 
consistency depends on the particular wave function and the r.m.s. transverse momentum 
chosen. Self-consistency is defined such that 50% of the result are accumulated in regions 
where a 2 s is smaller than 0.5 (Sec. V). 

Comparing our theoretical results with the data, it turns out that they fall short by 
at least 50%. This is true not only for the COZ DA, (which we have exemplarily used to 
facilitate comparison with previous works) but actually for the whole spectrum of amplitudes 
determined in |30|j3l|1 and found to comply with the COZ sum-rule requirements. Depending 
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on the actual value of the r.m.s. transverse momentum, the reduction of the perturbative 
contribution may be even stronger than 50%. 

The fact that in all considered cases the self-consistently calculated perturbative contri- 
bution to the proton form factor fails to reproduce the existing data, is perhaps a signal 
that soft contributions (higher twists) not accounted for so far by the modified HSP should 
be included. Such contributions comprise, e.g., improved and/or more complicated wave 
functions, orbital angular momentum, higher Fock components, quark-quark correlations 
(diquarks), radiative corrections to the quark and gluon condensates, quark masses, etc. 
Also remainders of genuine soft contributions, like vector-meson-dominance terms or the 
overlap of the soft parts of the wave functions (Feynman contributions), may still be large at 
accessible momentum transfers. The rather large value of the Pauli form factor F 2 around 
10 GeV 2 , as found experimentally [|1^], indicates that sizeable higher-twist contributions 
still exist in that region of momentum transfer ||13|| . One may suspect similar or even larger 
higher-twist contributions to the helicity non-flip current matrix element controlling F\ and 
Gm- Large (perturbative) higher-order corrections to the hard-scattering amplitude cannot 
be excluded as well, since their size has not yet been estimated. In analogy to the Drell-Yan 
process, these corrections may be condensed in a K-factor multiplying the leading-order 
perturbative result. However, with our choice of the renormalization scale, the K-factor is 
expected to be close to unity. At least for the case of the pion form factor, calculations of the 
K-factor to one-loop order exist |T6|JT7|1 , which indicate that choosing the renormalization 
scale analogously to ours, the value of the K-factor is indeed close to unity. 

In conclusion we note that it was not our primary aim to use the modified HSP to obtain 
best agreement with the data, although from our point of view this scheme represents a 
decisive step towards a deeper understanding of the electromagnetic form factors. In the 
present work the focus has been placed on theoretical problems, overlooked previously. 
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TABLES 



Table 1. Expansion coefficients for selected DAs, taken from [8] and [3C.31]. Our notation is 
adopted from J?]]. The {7™} are related to the anomalous dimensions of trilinear quark operators. 
The associated r.m.s transverse momentum and oscillator parameter for each model wave function, 
normalized via P^ q = 1, are shown. 



n 


7n 


£ n (COZ u P) 


£?n(COZ) 


^(COZ ?*) 


5n(het) 


-B n (low) 


1 


20/9 


3.2185 


3.6750 


3.5268 
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FIGURES 

FIG. 1. Illustration of gluonic radiative corrections to the proton magnetic form factor in the 
axial gauge. 



FIG. 2. The exponential of the Sudakov function s(£j, bi, Q) vs. and &/Aqcd for Q = 30Aqcd- 
In the hatched area the Sudakov function is set equal to zero according to Li's requirement f23|. 



FIG. 3. The Sudakov factor e _Sl vs. &iAqcd and ^Aqcd evaluated at Q = 30Aqcd, and 
x\ = Xi = 0.9, X2 = X3 = x' 2 = £3 = 0.05 assuming a linear quark configuration (b\ and 62 
are parallel to each other). The upper and lower figures correspond to the "L" and "MAX" 
prescriptions, respectively. 

FIG. 4. The proton magnetic form factor as a function of 6 c Aqcd- The curves shown are 
obtained at Q = 30Aqcd for the COZ DA. The solid line corresponds to the "MAX" prescription 
including evolution. The dotted (dashed, dashed-dotted) line represents results using the "H-SS" 
("MAX", "L") prescription ignoring intrinsic k± and evolution. 

FIG. 5. The influence of the intrinsic transverse momentum on the proton magnetic form factor. 
The curves shown are obtained for the COZ DA by imposing the "MAX" prescription including 
evolution. The solid line represents the results without &j_-dependence, whereas the dashed and 
dotted lines are obtained with < k\_ > l l 2 = 271 MeV and 600 MeV, respectively. 



FIG. 6. The proton magnetic form factor vs. Q . Data are taken from [15|. The Gr M data are 



represented by black circles, whereas those for Ff are indicated by open circles. The theoretical 
results are obtained using the "MAX" prescription including evolution and normalizing the wave 
functions to unity. The shadowed strip indicates the range of predictions derived from the set of 



DAs determined in [3C,31] in the context of QCD sum rules (see text). The solid (dashed, dotted) 



line corresponds to the COZ (heterotic, optimized COZ) DA. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405340vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405340vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405340vl 
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This figure "figl-6.png" is available in "png" format from: 
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